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By carrying out Monte Carlo simulations, we study step bunching during solution growth. For 
simplicity, we consider a square lattice, which represents a diffusion field in a solution, and express 
the diffusion of atoms as the hopping of atoms on the lattice sites. In our model, we neglect the 
fluctuation along steps. An array of steps is expressed as dots on a one-dimensional vicinal face. 
Step bunching occurs in the case of step-down flow. In previous study (M. Sato: J. Phys. Soc. Jpn. 
79 (2010) 064606), we studied step bunching with a slow flow and showed that the width of the 
fluctuation of step distance increases as t a with a = 1/3 in the initial stage. In this paper, we carry 
out simulations with a faster flow. With a faster flow, the width of the fluctuation of step distance 
first increases as t a with a — 1/2, which is larger than the exponent with a slow flow, and then an 
interval during which the width increases as t 1 ^ 3 appears. 



PACS numbers: 



I. INTRODUCTION 



On a vicinal face, which consists of an equidistant train 
of straight steps, the steps show two types of instability. 
One is step bunching, which is the instability of step dis- 
tance, and the other is step wandering, which is the insta- 
bility along steps. These instabilities are often observed 
during growth and melting. 

During solution growth, step wandering is caused by 
the flow in a solution. According to a linear stability 
analysis [l[ , step wandering is induced by a step-up flow, 
in agreement with experimental results |2| . The flow in a 
solution also causes step bunching. Chernov and cowork- 
ers [H-IH assumed that the step density is sufficiently high 
that the vicinal face can be treated as a linear sink of 
atoms, and they studied the step bunching caused by the 
flow in a solution. Linear stability analyses and nu- 
merical studies [H-IH have shown that a growing vicinal 
face is unstable in the case of a long-wavclcngth fluctua- 
tion owing to a step-down flow. The time evolution of an 
unstable vicinal face was studied by Brcdikhin and Mal- 
shakova Q . According to numerical simulations based on 
a nonlinear equation, a quasi-regular array of high step 
density, which is similar to that of bunches observed in 
experiments [|| 0] , is formed on an unstabilized vicinal 
face. 

In previous studies the modulation of step den- 

sity during step bunching was investigated but the mo- 
tion of each step has not been sufficiently studied. Re- 
cently, wc made a simple discrete model and, by carrying 
out a Monte Carlo simulation Q, studied the motion of 
each step during bunching. In the model, we assumed 
that the steps are straight and we expressed the steps as 
points on a one-dimensional vicinal face. To represent 
the diffusion field in a solution, we introduced a square 
lattice and treated the diffusion in a solution as the hop- 
ping of atoms on the lattice. In the paper 0] , we carried 
out a Monte Carlo simulation and showed that the vicinal 
face is unstable in the case of a long-wavelength fluctu- 
ation. The width of fluctuation of the step distance, w, 



increases as t a with a = 1/3. In the simulation, we used 
low flow rates. If the flow rate is changed, the time evolu- 
tion of the fluctuation of step distance and the behavior 
of steps during bunching may be changed. 

In this paper, we carry out Monte Carlo simulations 
with a flow in a solution whose flow rate is higher than 
that in previous paper @ and study how the time evo- 
lution of bunching changes. In ClTI we introduce our 
simplified model. In £ )III[ we report the results of Monte 
Carlo simulations. The width of the fluctuation of step 
distance increases with time in the initial stage and then 
becomes saturated. The growth law of the width of the 
fluctuation and the dependence of the saturated width 
on supersaturation are studied. In mV[ we summarize 
the results and give a brief discussion. 



II. MODEL 

Using a very simplified model (Fig. [1]) we carry out 
Monte Carlo simulations. In the model, a square lattice 
in a solution is considered. The diffusion of atoms in a 
solution is expressed as hopping on the lattice sites. The 
height of steps is neglected, a vicinal face is expressed as 
a line, and steps are treated as points on the line. The 
x-direction is the step-down direction and z — is the 
altitude of the vicinal face. 

In our model, the active atoms are solid atoms forming 
steps and atoms in solution. In each Monte Carlo trial, 
we randomly choose one active atom. When an atom 
in a solution is chosen, a diffusion trial is carried out. 
The flow in a solution is expressed as the asymmetry of 
the diffusion probability [§, |T(| ■ When a flow is prarallcl 
to the x-direction, an atom at site moves to site 

(i, j ± 1) with probability 1/4 and to the site (i ± l,j) 
with probability, p± = (l±pfl ow )/4. The parameter paow 
represents the effect of the flow in the solution. With in- 
creasing parameter paow, the effect of the flow increases. 
The time is increased after each diffusion trial. In order 
that the diffusion coefficient D is equal to 1, the incre- 
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FIG. 1: Schematic figures of (a) a vicinal face and (b) our 
model. Squares and circles represent solid atoms and atoms 
in solution, respectively. Dark squares are the solid atoms 
that form steps. Arrows show the positions of steps. 



merit of the time after each diffusion trial is given by 
At = l/4iV s , where N s is the number of atoms in the 
solution. 

To take account of the viscosity in the solution, we as- 
sume that pflow is proportional to the altitude from the 
vicinal face. When the z-coordinate is given by z = j, 
p fl ow is expressed as p flow = Pf\ ovr (j - !)• Tnu s, near 
the vicinal face, the effect of the flow in the diffusion 
trial vanishes. The sign of Pf\ ow shows the direction of 
the flow. When Pf\ OVJ is positive, the flow is in the step- 
down direction, and the flow is in the step-up direction 
for negative pnaw- If the altitude is sufficiently high that 
Mow > Inflow < —1) j we set the hopping probabilities 
p + and p_ to be 1(0) and 0(1), respectively. In general, 
we should consider a flow obeys that boundary-layer the- 
ory. However, according to the Blasius's solution, the ve- 
locity of a smooth flow near a flat surface is proportional 
to the distance from the surface. Thus, if the effect of 
the flow far from the surface is neglected, our model is 
valid as a simplified model. 
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FIG. 2: Schematic figures of solidification process in (a) a 
vicinal face and (b) our model. 

The solidification probability p+ and the melting prob- 
ability p_ are given by [ll[ 



1 



P± 



(1) 



where <j) is the chemical potential gain upon solidification. 
A solidification trial is carried out after a diffusion trial. 
If the position of the atom is (if , 1) and that of the step 
is x = if — 1, the step advances to x = if and the atom 
in the solution vanishes after solidification (Fig. 0. A 
melting trial is carried out when the chosen active atom 
is an atom forming a step. If the melting of a step atom 
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FIG. 3: Schematic figures of melting process in (a) a vicinal 
face and (b) our model. 



occurs at a step edge x = i- u the step recedes to x = i\ — 1 
and an atom is expelled to site («i,l) (Fig. [3J. In our 
algorithm, only one atom can occupy a site. Then, if site 
(ii, 1) is already occupied by an atom, the melting trial 
is not carried out (Fig.@|. 
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FIG. 4: Forbidden melting process in our algorithm. Only 
one atom can occupy one site. 

We assume that the kinetic coefficient at steps is suffi- 
ciently large that the detailed balance stands up at step 
positions. Namely, the concentration of atom is in equi- 
librium at steps. The equilibrium concentration in a so- 
lution, c cq , satisfies 



CcqPs 



(1 - C eq )p n 



(2) 



From eq. @, the equilibrium adatom density is given 
by [Hi 



1 
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(3) 



In our model, when the altitude is larger than the crit- 
ical value z max , the concentration of atoms becomes a 
constant value c^. If is higher than c° q , solidification 
frequently occurs and the steps advance. In the opposite 
case, melting occurs more frequently than solidification 
and the steps recede. 



III. RESULTS OF SIMULATION 

We previously studied the motion of steps during step 
bunching using the same model Q. In the previous pa- 
per, we showed that a step-up flow stabilizes the vicinal 
face and that the step bunching occurs with a step-down 
flow. Thus, hereafter, we investigate the time evolution 
of steps with a step-down flow. In the previous study [|| , 
we studied the motion of steps with low flow rates. In 
an early stage of simulation, a vicinal face seems to be 
unstable for fluctuations with a long wavelength. The 
step density gradually modulates and small bunches are 
formed. When step bunching is caused by the drift of 
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adatoms [l2( , the formation of bunches changes with the 
velocity of the drift flow. A similar change may be gen- 
erated in the present model. In the previous study ||, 
the flow rate is sufficiently slow that the width of the re- 
gion with a shear flow is 10 — 100 lattices. In the present 
case, we consider that the shear rate is sufficiently large 
that the width of the region is narrower than 10. We 
carry out simulations with a fast flow and study how the 
growth law changes. First, we study the time evolution 
of the width of the fluctuation of step distance, w, which 
is defined as 



JV 



(4) 
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N is the number of samples, n s is the number of steps, 
is the width of the ith terrace in the jth sample, 
and I is the average terrace width. 
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FIG. 6: Initial stage of time evolution of width of step dis- 
tance with (a) Coo = 0.2 and (b) Coo = 0.4. 




FIG. 7: Initial stage of time evolution of width of step dis- 
tance with pflow = 0.005. The other parameters are the same 
as those in Fig. Ha). 



FIG. 5: Typical time evolution of width of fluctuation of step 
distance. The parameters are 0//cbT = 3.0, Pfl ow = 0.3, and 
Coo = 0.2. 

Figure [5] shows a typical time evolution of the width 
of the fluctuation of step distance. We used the param- 
eters 4>/kBT = 3.0, pflow = 0.3, and Coo = 0.2. The 
system sizes in the ^-direction and ^-direction are given 
by L x = 512 and L z = 128, respectively. Since n s = 16, 
the average step distance is I = 32. The data is averaged 
over 100 samples. In the initial stage, w increases with 
time. Then, the saturation of w occurs. Although the 
time evolution seems to be similar to that with a slow 
flow, we study the time evolution in detail. 

Figure [5] shows the initial stage of the time evolution. 
The density of atoms far from the vicinal face is given 
by Coo = 0.2 in Fig.^a) and Coo = 0.4 in Fig.^b). The 
other parameters are the same as those in Fig. [5J As 
shown in Fig. [51 the exponent changes with two stages. 
In the first stage, the exponent is given by a = 1/2, which 
is larger than that with a slow flow. In the second stage, 
the growth of w becomes slow and the exponent is given 
by a = 1/3, which is the same as that with a slow flow. 
In previous study [§[, we studied the time evolution for 
various slow flows. In Fig. we show the time evolution 
with Pfl ow = 0.005 as a sample of the time evolution with 
a slow flow. In an early stage (10 2 < t < 10 4 ), no change 
in the exponent occurs and a = 1/3. 



The time at which the exponent a changes from 1/2 
to 1/3 is about t w 10 3 in Fig. Efa) and t « 3 x 10 3 in 
Fig.[6][b). The time interval with a = 1/2 in Fig.[6][a) is 
shorter than that in Fig.^b). This difference is probably 
caused by the increase in the supersaturation. However, 
the time at which the saturation of width occurs hardly 
changes irrespective of the supersaturation. Thus, the 
decrease of time interval with a = 1/3 is probably caused 
by the decrease in the supersaturation. If we use a larger 

, the time interval in which the exponent is a = 1/3 
may vanish. 

Figure [8] shows the time evolution of step positions. 
The parameters are the same as those in Fig. [SI The 
value of a is related to the motion of steps. In the initial 
stage, in which w increases as t 1 / 2 , single steps gather 
and form small bunches. The vicinal face separates into 
bunches and terraces. Steps do not exist on the terraces. 
With increasing supersaturation, the single steps advance 
rapidly and tight bunches are formed. In the later stage, 
in which w increases as t 1 / 3 , both the growth of bunches 
by the coalescence of small bunches and the separation 
of steps repeatedly occur. 

In the final stage, in which the saturation of w oc- 
curs, the collision of bunches, the temporary breaking of 
bunches into steps and small bunches, and their recoales- 
cence repeatedly occur Q . To determine the dependence 
of the saturated width on the supersaturation, we 
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FIG. 8: Time evolution of step positions with (a) Coo = 0.2 
and (b) c x = 0.4. The parameters are the same as those in 
Fig.® 



carry out simulations with various . Figure |H] shows 
the dependence of on the supersaturation {crx, — c oq ). 
Since we use (jj/k^T = 3.0, the equilibrium adatom den- 
sity is given by c eq = 4.7 x 1CP 2 . We use the flow rate 



Pflow = 0.1. With increasing supersaturation, 
creases. The dependence is given by Woo ~ (coo — c eqJ 
The value of 7 is approximately given by 7 = 3/20, which 
is similar to the dependence of Woo on pnow 
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When the flow is rapid, the width of the fluctuation of 
step distance increases with time as t a and a changes 
with two stages. In the first stage, a is equal to 1/2, 
which is larger than that with a slow flow. In the second 
stage, a is 1/3, which is the same as that with a slow 
flow. The time interval with a = 1/2 becomes longer 
with increasing supersaturation. 

The change in the exponent is related to the motion 
of the steps. When a = 1/2, bunches grow by the co- 
alescence of single steps. A vicinal face is separated to 
form a wide terrace and small bunches consisting a few 
steps. The exponent a = 1/2 means that the correlation 
between steps is weak and that the steps move randomly. 
The situation is similar to the roughing of an interface 
by ballistic deposition [l3j]. When a = 1/3, the step 
separation of steps from bunches repeatedly occurs. Sin- 
gle steps are present on large terraces. In a later stage, 
the separation of steps is more frequent than that in an 
early stage. Owing to the single steps on large terraces, 
the average terrace width decreases. Thus, the growth 
of the width of fluctuation becomes graduall and a is 
lowered. In the last stage, the saturation of the width 
w occurs. The saturated value increases as (Ac) 7 
with 7 3/20. The dependence of on the supersat- 
uration is the same as that of the flow rate Q. We are 
currently attempting to clarify how the exponents can be 
determined. 

In this paper, we assumed that the steps are straight. 
Thus, step wandering, which is the instability due to the 
fluctuation along steps, cannot be treated in our model. 
To study step wandering, we intend to carry out simula- 
tions using a two-dimensional step model. 



FIG. 9: Dependence of Woo on the supersaturation, where 



Ac = (co 



The flow rate is given by pfl ow = 0.1. 
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